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I. INTRODUCTION 



The complexation of oppositely charged polyelectrolytes in different ionic conditions is 
an interesting problem of biological relevance. Moreover, the association of two oppositely 
charged linear chains has implications in the design of new materials with unique proper- 
ties such as multilayer polyelectrolytes^'^, and carrier gels for drug delivery^. Industrial 
applications include uses as coatings, flocculants and absorbents. ^'^ 

Charged chains can associate via many mechanisms given the large number of length 
scales involved in these mixtures. The properties of the complexes are strongly dependent 
on physical parameters such as salt and monomer concentrations, and on the chemical config- 
urations of the chains such as their charge densities, persistence lengths, degrees of polymer- 
ization, and nature of the interactions between the charge groups along the backbone. 

Complexes formed by chains with low charge densities have been analyzed using lin- 
earized models. ^^^^^ These models are applicable to describe solutions with strongly hydrated 
charge groups of positive and negative charges, which interact weakly with each other. In 
these studies ion condensation effects^°'^^ can be neglected if the hnear charge density is suffi- 
ciently low. The complexation in these systems occurs via collective charge fluctuations. 
With addition of salt the complexes dissolve due to screening of electrostatic interactions. 

In certain polymer mixtures, even when the charge density is low, non-linear effects 
are important if the charge groups are strongly interacting. In these mixtures the oppo- 
sitely charged groups, when placed at short separation distances, can be locally dehydrated, 
and act as localized short-range crosslinks between oppositely charged groups along the 
chains. The formation of these links generates a thermoreversibly crosslinked solution^^"^^ 
of oppositely charged chains. The number of crosslinks formed in equilibrium in such ther- 
moreversibly associating chains is given by the law of mass action with an effective asso- 
ciation constant. In uncharged systems the association constant depends exponentially on 
the strength of the short range attraction between the reactive groups. For charged reactive 
groups, however, there is an additional term in the association constant due to the electro- 
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static contribution to the free energy resulting from the collective charge fluctuations. Since 
association occurs between the ions belonging to oppositely charged chains, this additional 
electrostatic term always increases the rate of association. The importance of non-linear 
association of charges has been recently recognized in polyelectrolyte adsorption^^'^^ and 
multilayer formation^^. 

In this work the electrostatic interactions are described using a two-fold approach. 
The strongly non-linear short-range interactions between oppositely charged groups are ac- 
counted for by including strong correlations between the chains. The long-range electrostatic 
interactions (which are weak for weakly charged polyelectrolytes) are accounted for in a lin- 
earized way by computing the fluctuations of this correlated solution of crosslinked charged 
chains using a generalized Debye-Hiickel approach (Random Phase Approximation). ^^'^^ The 
electrostatic free energy depends on the structure of charged polymer clusters, while the 
cluster distribution depends (through the modified law of mass action) on the electrostatic 
free energy. Therefore, we determine here the number of formed crosslinks by evaluating 
self-consistently the electrostatic contribution from the collective charge fluctuations of a 
crosslinked system of charged chains. 

The degree of hydrophobicity of the chain backbone modifies the thermodynamics of the 
solution. We consider non-selective solvents, where the degree of compatibility is the same 
for both the positively and negatively charged chains. We investigate how the degree of 
hydrophobicity infiuences the properties of complexes formed by electrostatic interactions. 

The paper is organized as follows. In Section 11 we describe the model and approximations 
used, and derive the free energy of the solution; the details of the derivation of the correlation 
function of the crosslinked system required to determine the electrostatic contribution to 
the free energy is given in Appendix A. In Section III we discuss how different system 
parameters influence the properties of the formed complex and its response to addition of 
salt. The conclusions are given in Section IV. 



3 



II. THE FREE ENERGY OF THE SEMI-DILUTE SOLUTION 

A. Model 

In this section we calculate the free energy of a homogeneous semidilute solution of op- 
positely charged polyelectrolytes. For simplicity in this work we consider only the absolutely 
symmetric case. Positively and negatively charged chains are present in the solution in equal 
concentrations, and have the same physical properties except the sign of the charge. The 
chains have degree of polymerization N = Ni = N2. The fraction of charged monomers 
on both types of chains is equal to /. We consider only weakly charged chains for which / 
is small enough, so that electrostatic energy of adjacent along the chain charges is smaller 
than the thermal energy. The number concentrations of positively and negatively charged 
monomeric units are pi — P2, so that the total concentration of monomeric units in the so- 
lution is p = 2pi//. Each charged monomeric unit releases a monovalent counterion. A 1:1 
low molecular salt can also be present in the solution with concentration of positively and 
negatively charged ions given by Ps+ and Ps-, respectively, and the total concentration of 
salt ions ps — '2ps+ — '^Ps-- Since the effect of counterions is equivalent to addition of salt 
we include the counterions in the salt. 

We describe strong electrostatic interactions between oppositely charged monomeric units 
as thermoreversible bond (crosslink) formation. This strong non-linear interaction should 
be treated differently from the long-range electrostatic part which is treated within the 
Random Phase Approximation (RPA).^^'^^ Another reason why formation of ionic pairs 
should be considered in addition to Coulomb interactions is because formation of ionic pairs 
can proceed with the rearrangement of the solvation shell of charges on polymer. In this case 
the effective dielectric constant between charges in a pair can differ significantly from the bulk 
solvent dielectric constant. Thus, the bond energy in a pair can be quite high. A natural 
model to describe such ion-pairing is reversible association between charges of opposite 
signs. For simplicity we assume that only pairs can be formed, with the absolute value of 
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the reduced bond energy e — \E\ /kT {k is the Boltzmann constant, T the thermodynamic 
temperature), which gives rise to association constant cu = e^. Note that bond formation is 
different from short-range van der Waals attraction in that it has the saturation property, 
that is, once a bond between two given ions is formed they do not interact with any other 
ions. 

We write the free energy of the solution of associating polyelectrolytes in the following 
form 

F = Fref + Frpa (1) 

where the first term is the free energy of the reference neutral system (but with short-range 
interactions) and the second term is the contribution of electrostatics. The electrostatic part 
Frpa is calculated within the RPA, which is a linear theory equivalent to the Debye-Hiickel 
approximation (i.e., to the linearized Poisson- Boltzmann equation). For the reference free 
energy we use the Flory-Huggins mean-field approximation 



re/ 



Fid + Ffh (2) 



which includes the ideal entropic and enthalpic terms. In our associating system the ideal 
term is the free energy of ideal gas of all possible clusters {C} with appropriate statistical 
weights oj{C) 

Fid /^x, P{C) Ps- Ps+ /„s 

= P{Q In — + ps~ In h ps+ In (3) 



kTV f^^'^'^ eu;{C) ' e ' e 

as well the entropy of the ideal gas of salt ions (counterions are also included here) . As has 
been shown in the refs 22, 23 the equilibrium concentrations p{C) can be obtained using a 

diagrammatic technique and the free energy of associating chains can be written as 



^ P^^^ ^ = ^ In p + p/ [(1 - r) ln(l - r) + r In r] - ^ In 
{C} eLU[C) 1\ Z 



2e 



(4) 



Here conversion F is the fraction of polymeric ions in pairs (see Appendix A), which is to 
be found from subsequent minimization of the total free energy of the solution. The volume 



of the monomeric unit v (which we for simphcity assume to be equal to v — b^) is used 
to approximate the internal partition function of the crosslink Z cross — ve^. Alternative 
combinatorial derivation of (4) can be performed along the lines of ref 32. 

We assume that the polycation and polyanion backbones have identical short range 
interaction with solvent. The interaction free energy in (2) is assumed to be given by the 
Flory-Huggins form 

^|^ = (l-0-0,)ln(l-0-0,)+x0(l-0) (5) 

where V is the volume of the system, = pb^ is the total polymer volume fraction, and 
4>s = b^ps = '2b^Ps- is the total volume fraction of salt ions (and also counterions) . The 
first term in (5) stems from hardcore repulsion; the second one from short-range attraction, 
whose strength is characterized by the parameter x- We will analyze both the cases of good 
and marginal to bad solvent. Note that, in contrast to previous works, ^^'^^ interactions of 
backbones with solvent favor complexation under bad solvent conditions. 

Adding up the two contributions, the free energy Fref of the reference neutral system 
reads 



^ 1 ln0 + 0,ln0, + 0/ [(1 - r) ln(l - r) + rinr] - ^ In 
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+ 



+(1 - - 0,) ln(l -<l>-<f>s)+ X0(1 - (/>) (6) 



B. Electrostatic free energy: Random Phase Approximation 

Due to electroneutrality the electrostatic contribution F^pa in the total free energy (1) 
is due to fluctuations of charge concentration, which is calculated within the Random Phase 
Approximation (see ref 23 for details): 
Frpa V f cPq 



I 



In (det (I + G{q)V{q))) - Y^PiUM 



(7) 



kT 2 J (27r)3 

where I = is the unitary matrix, G(g) is the correlation function matrix of the reference 
neutral system and U(g) is the matrix of Coulomb interactions. The sum runs over all 
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charged components of the system (co-ions and salt ions). The last term in (7) is the self- 
energy of pointlike charges. The correlation function matrix G{q) can be in turn obtained 
within the RPA as^''"^^ 



(8) 



where g{q) is the structure correlation matrix. It characterizes correlations of density due to 
existence of different clusters in the system, but does not include interactions. The matrix 
g(g) for our symmetric system has the form 



9ij 



(9) 



gu 912 

5-12 gu 

ps- 
\ p,+ J 

The polymeric correlation function gn and gi2 are calculated in Appendix A using a dia- 
grammatic approach. The interaction matrix c{q) describes short-range interactions (free 
energy Fpn) and its components are given by 

1 



SiSj - 2xPiPj 



(10) 



where we introduced the volume fraction of solvent $ = 1 — — 0s and two auxiliary vectors 



Si = {1,1,1,1} 
p, = {l,l,0,0} 



(11) 
(12) 



Using the vector of valencies Ci, the Coulomb interaction matrix can be written in the 
following form 



Uij{q) = eiCjUiq) 

e, = {1,-1, 1,-1} 



(13) 
(14) 



which allows us to simplify the expression under the logarithm in (7) 



det (I + G{q)V{q)) = 1 + U{q) ^ G,,(?)e,e,- (15) 

Using formulae (10) and (12) for the correlation functions gij and Qj we now can obtain Gij 
in accordance with (9). Substituting the result into (15) we finally obtain 

det (I + G{q)V{q)) = 1 + U{q) {2gn - 2gu + 2p,_} (16) 

It is remarkable that this result is the same as if we had neglected the short-range interactions 
(the matrix Qj), that is, if we had put G^j = g^j in determinant (15). Note, however, that 
this simple result holds only for our case of a symmetric system (described by matrix g^j) 
and for symmetric long-range (matrix U(g)) and short-range (matrix c) interactions. 

The structure correlation functions g^j are calculated in Appendix A (see (A19-A20)). 
We reproduce them here for convenience 

/ ^l + (^')'M?) 

..(,)^P..(«)j^j^ (18) 

The functions g{q) and h{q) are defined by (A6) and (A8) in Appendix A. The effective 
conversion V is defined as Te~'^^^^/^ in (A18), with the bare conversion F defined as the 
fraction of charged monomers participating in crosslinks (see eq A14). It is important 
to note that the correlation functions (17-18) are calculated for an ideal thermoreversibly 
associating system, in which no other interactions except crosslinking are present (in our 
case no electrostatic and hydrophobic interactions).^^ 

Substituting these expressions into (16) we can rewrite the free energy (7) as 



Frpa V f d^q 



(19) 



kT 2 J (27 
where pf = 2pi = 2p2, and Ps = 2ps_ = 2ps+. 

To be able to evaluate Frpa we need to specify a suitable form of the interaction potential 
U{q). For bare Coulomb interaction we have 

kT ekTr r ^ ' 

UM_ fd^re^-^Uc{r)^^4 (21) 



kT J ' ?2 



Here we introduced the Bjerrum length I — ql/{ekT), where Qe is the electron charge and 
e the dielectric constant of the solvent. In order to take into account the influence of the 
hardcore of the ions on the electrostatic contribution Frpa we use a modified Coulomb 
potential given by 

^^-fl-e-^/") (22) 

kT q^{l + q%^) ^ ' 

where the bond length h is for simplicity taken to be the size of the ions. At large distances 
(r ^ h) the modified potential becomes the pure Coulomb potential (20-21). However, 
at r = the modified potential attains a finite value, while the original Coulomb poten- 
tial diverges. Thus we phenomenologically include the impenetrability of the ions within 
the RPA formalism, which is originally formulated for pointlike ions. The RPA with the 
modified potential (22-23) has been shown to successfully describe the phase diagrams of 
polyelectrolytes^^ and of the low-molecular system of charged dumbbells.^^ Furthermore, 
this potential has been successfully used in the liquid state approaches. 

Substituting the modified potential (23) into (19) we obtain the final expression for the 
electrostatic free energy 



kT 2 J (27r 



(24) 



Let us make several comments on the structure of F^pa- The polymer structure correlation 
functions (17-23) diverge at the gelation condition r{Nf — 1) — 1, which is simply due to the 
fact that gelation corresponds to the formation of an infinite cluster. It is remarkable that 
the electrostatic free energy (24) has no corresponding singularity at the gelation structural 
transition. This is due to the charge symmetry of the considered system. Indeed since in 
our case association is possible only between oppositely charged chains (which carry the 
same amount of charge) the infinite cluster is by construction neutral, therefore does not 
contribute to Frpa- (It can be shown that for any asymmetric system (asymmetry of N , p 
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or /) the infinite cluster is charged and, accordingly, expression for Frpa has a singularity 
at and beyond the gelation transition. However, this unphysical singularity is an artifact of 
our simplified description of gel structure.) 

Let us look at the limiting cases of eq 24. By putting F' = we regain the well-known 
expression for free unassociated chains, and if we put F' = 1, the polyelectrolyte chains do 
not contribute to Frpa- However, since F' is only the effective conversion: F' = Ye~'^^^^^^ 
the equality F' = 1 is possible only when F = 1 and q — Q. On all finite lengthscales {q ^ 0) 
charges in crosshnks still contribute to the electrostatic free energy Frpa- This is natural, 
since in our model of crosslinking the two opposite charges do not annihilate, rather they are 
considered as separate charges with Gaussian correlations between them, which leads to the 
emergence of effective conversion F', instead of bare conversion F. Because the charges do 
not annihilate when the ionic pairs are formed, we subtract the self-energy of all ions present 
in the system (last term in (24)), regardless of whether they are free or form crosshnks. 

The chain correlation functions g{q) and h{q) are defined by (A6) and (A8). The function 
g{q) can be easily calculated in the continuous limit, the result being the well-known Debye 
structure function. However, we also need the correct limit of pointlike ions at 5 — > oo in 
(24), since we subtract the self-energy of all ions Therefore, we choose a simple interpolation 
form for the chain structure function g{q) 

Hi) = 9(9) - 1 (26) 

which gives the correct limit at g = 0, has the scaling of a Gaussian chain at N~^/'^ <^ qb <^ 
/^/^, and reproduces pointlike ions at qb /^/^. 

C. Minimization of the free energy 

The total free energy (1) is given by the sum of F^e/ in (6) and Frpa in (24). However, 
this is only the virtual free energy of a system with a given value of conversion F. To obtain 
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the equilibrium free energy F(req) we need to obtain the equihbrium conversion Feq, by the 
minimization of F(r): 

^^ = (27, 
Using (6) and (24) we obtain the following equation for F 

^ ^ exp [e + pLRPA{r)] (28) 



(1 - r)' 2 

This equation has the general structure of the law of mass action. A noteworthy feature of 
(28), however, is that the energy gained from formation of a crosslink consists of the bonding 
energy e and the energy gain /irpa resulting from the long-range electrostatic attraction of 
the polymer chains described by Frpa'- 

AttJ 1 

Of course the energy fiRPA depends on the thermodynamic state of the solution, and thus on 
r, therefore equations (28) and (29) are to be solved simultaneously to obtain T{(f)). Since 
in our model crosslinking takes place only between oppositely charged chains /iRPA(r) is 
always positive, that is, the electrostatic attraction, along with the specific binding energy 
e, always promotes crosslinking. From the numerical solution one obtains that /XijpA(r) is a 
monotonically decreasing function of F for all (f), which is explained by the fact that as more 
ions associate they contribute less to the long-range attraction (which can be seen directly 
from Frpa in (24)). 



III. RESULTS AND DISCUSSION 



Due to its complete symmetry, only macroscopic phase separation is possible in the 
considered system. In our ternary incompressible system of polymer, salt and solvent we have 
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two independent components, which we choose to be polymer and salt. When macroscopic 
phase separation (precipitation) occurs the two coexisting phases differ in concentrations of 
both polymer and salt. Thus, generally speaking, we have to calculate the phase diagrams 
of a ternary incompressible system. However, we are mostly interested in two aspects of 
the precipitation process: the influence of different competing complexation mechanisms on 
the density of the formed precipitate (studied in the next section) and determination of the 
conditions of solubihty of the complexes with addition of salt (section IIIB). In both of 
these cases we can make speciflc additional assumptions which allow us to investigate the 
mentioned problems in a simple and clear manner. 



A. Density of precipitate: effect of crosslinking and van der Waals attraction 

In our model we consider three complexation driving forces: long-range electrostatic 
attraction between co-ions, strongly non-linear short-range attraction leading to ion- 
crosslinking and van der Waals attraction between all monomeric units. In this section, 
we look at the density of the formed complex 0, in particular how depends on the relative 
importance of the three complexation factors as well as how the density is influenced by 
the addition of salt. We can signiflcantly simplify the analysis if we make the following 
two assumptions (similar to the ones employed in our previous work).^^ First, we assume 
an inflnite degree of polymerization N = oo. Second, the total concentration of polymer 
chains in the whole solution (system) is assumed to be small. Since the entropy of the chains 
represents the only driving force for dissolution of the polymer chains from the precipitate, 
the flrst assumption amounts to assuming zero polymer concentration in the supernatant 
(which for a flnite N would be a polymer-poor phase) . The second assumption is equivalent 
to assuming that the salt volume fraction in the supernatant is equal to the salt volume 
fraction in the whole system, which we thus denote simply as 4>s- Note that the salt volume 
fraction in the precipitate (j)^^ can differ considerably from 0^, which, as we show below, has 
a signiflcant effect on 0. Given our assumptions, and (f)^^ can be found by equating the 
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pressure and the chemical potential of the salt in the coexisting phases: 

p(0 = O,0,)=p(0,0(^')) 
/x,((/. = O,0,)=/x,(0,0(f)) 

For convenience here and in the following we use a dimensionless equilibrium free energy 
density defined by ^ = F{Vf,q) /{kTV /v). The equilibrium free energy F(Teq) is obtained 
from the minimization (described in section II C) of the total free energy, which according 
to (1) is given by the sum of (6) and (24). We solve equations (32-33) numerically to obtain 
and considering 0^ as a parameter. The results are given in Figures 1-3. 

In Figure 1 we assume no van der Waals attraction x = and plot results for varying 
bonding energy e — E/kT. (Experiments on polyelectrolyte adsorption^^'^^ and multilayer 
formation^^ are consistent with the value of binding energy e varying between £ = 3 and 
e = 7.) We set the Bjcrrum length / = 3, and the fraction of charged monomers / = 0.1. 
The dependence of on Z and / for the case of complexation without crosslinking was 
investigated in our previous work.^^ Similar to the results for that case, increases with 
increasing / and/or / for all 0s, e or x- With good precision one can obtain results for other 
values of / simply by linearly scaling with /. 

Figure 1(a) shows the polymer volume fraction in the precipitate as a function of the 
salt volume fraction in the system 0^. We see that for all values of s the complex density 
monotonically decreases with increasing salt concentration. This is of course due to Debye- 
Hiickel screening by salt, which makes the electrostatic attraction weaker (the term Frpa in 
the total free energy). Note that as Frpa becomes weaker the electrostatic binding energy /x, 
given by (29), also becomes smaller, thus both the long-range and short-range electrostatics 
become less efficient in forming the complex. This is illustrated in Figure 1(b), where we 
plot conversion F for curves of plot Figure 1(a). We see that conversion also monotonically 
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(32) 
(33) 
(34) 

(35) 



drops when more salt is added to the system, the behavior being similar for all values of 
£. With the increase of £, conversion F increases for all salt concentrations. However this 
does not directly translate into a denser complex as we can see from Figure 1(a). Indeed for 
0s ~ we see that depends non-monotonically on e , the feature which will be explored in 
detail in Figure 3. For large salt concentrations the density is larger for greater £, which is 
obviously due to increased crosslinking of chains. The presence of crosslinks manifests itself 
in a feeble dependence of on 0^ for larger values of e, which indicates indissolubility by 
salt of complexes formed primarily by specific short-range attractions. At the same time for 
small binding energies (such as £ = and £ = 3) in Figure 1(a) the density 4> drops rather 
abruptly and becomes very small, which indicates dissolution of the complex with addition 
of salt (this problem will be studied in the next section). The difference between the salt 
concentration in the precipitate (p^f^and that in the supernatant (pg (which , according to our 
assumption, is equal to the salt concentration in the whole system) is presented in Figure 
1(c). We see that for small e the precipitate is first enriched with salt, which (as we showed 
previously^*) is due to correlational Debye-Hiickel attraction. For small s with increasing 
salt concentrations becomes smaller than 0^ and then as the complex becomes very 
diluted for large cpg there is only a neghgible difference. For large e we have cf)^/^ < cpg for 
all salt concentrations. (Note that we plot the difference (p^^ — (pg in Figure 1(c), of course 
0^) increases with (pg) The depletion of salt was previously shown to be due to hardcore 
interactions.^* This depletion turns out to have a considerable effect on (p. In Figure 1(a) we 
plot with dashed lines the curves for e = 5 and e = 7 obtained from equation (32) with the 
assumption = (pg. The effect can be seen to be especially substantial for larger values 
of (pg. 

In Figure 2 we illustrate the effect of the x-parameter for the case when the effect of 
association is small {e — G). The complex is seen to become denser with increasing short- 
range attraction for all values of (pg. For x < 0.5 (good solvent) the volume fraction (p 
strongly decreases with increasing salt (which is an indication of dissolution of the complex 
for finite A^). However, for x > 0.5 (bad solvent condition) the density is seen to be negligibly 
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dependent on 0^ for large 0s, with the complex being stable with respect to addition of 
salt. Comparing Figures 1(a) and 2 we observe that behavior of with increasing (pg is 
qualitatively the same when the complex is formed by crosslinking (e) or van der Waals 
interactions (x). Note that we analyze here only marginally bad solvent conditions, for 
which we can disregard the possibility of necklace formation. ^'^ 

The relative strength of the long-range electrostatic attraction vs. crosslinking is demon- 
strated in Figure 3. We plot the precipitate density in the salt-free solution as a function 
of bond energy e. Different curves correspond to varying x, which spans good to marginal 
solvent conditions. Increasing e means that a greater number of charges form crosslinks. 
Indeed from the numerical solution we obtain that in all cases F monotonically increases 
with growing £. Thus increasing s physically means changing the driving force of complex- 
ation from long-range charge correlation to crosslinking attraction (numerically, at £ = 10 
the conversion F f5i 1). Interestingly, as we see from Figure 3, the density (although it 
exhibits a shallow minimum as a function of e) is rather insensitive to the value of £. Thus, 
ion-pairing and long-range correlations lead to polyelectrolyte complexes of similar density 
and the two mechanisms can be difficult to distinguish experimentally for salt-free systems. 
Curves for different x show qualitatively the same behavior, with (p increasing as the solvent 
worsens. 



B. Phase diagram: dissolution with addition of salt 

As can be seen from Figures. 1(a) and 2 for small e and/or x the density becomes 
very small with addition of enough salt, which is indicative of precipitate dissolution for a 
finite N. In the previous section we assumed A?" = oo, so the precipitate never dissolved. In 
this section we relax this assumption and look at the phase coexistence. 

Let us first investigate the effect of N on the phase coexistence. In order to simplify the 
presentation let us assume that the salt volume fractions in the polymer-rich and polymer- 
poor phases are the same. Thus we treat 0^ as a parameter and obtain the coexisting 
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polymer concentrations by constructing a common tangent to the equilibrium free energy 
density T . An example of resulting phase diagrams for varying N is shown in Figure 4 
(the values of all parameters are given in the plot). The area below the coexistence line 
(for a given corresponds to phase separation (complexation at low salt concentrations), 
above — to a homogeneous solution (the precipitate dissolves at high salt ^5). Owing to 
our simplifying assumption of equal 0s in the two phases, the tie lines giving the coexisting 
polymer concentrations are parallel to the 0-axis. We see that the dilute phase has neghgible 
polymer concentration even for rather small A'^, which justifies the assumption N — 00 oi 
the previous section. When enough salt is added to the solution the precipitate dissolves. As 
we can see from Figure 4 the salt concentration needed to dissolve the precipitate strongly 
depends on N and (as Figures 1 and 2 show) it also depends on e and x- 

In Figure 5 we determine the conditions on £, and A^ ensuring stability of the pre- 
cipitate at a given concentration of salt. For each curve (corresponding to a certain x) the 
precipitate exists at 0^ — 0.1 if the values of e and x lis in the area above the curve, and 
the precipitate is dissolved in the area below the curve. (The value 0^ = 0.1 is taken as 
an example and it applies to all curves). The increase of either of £, x, or A^ stabihzes the 
complex to addition of salt. We observe that the stabihty of the precipitate is quite sensitive 
to the values of parameters in the experimentally most relevant region £ pa 3, x ~ 0.5, and 
100 < N < 1000. The results of Figure 5 can be used for experimental design of complexes 
stable to salt. 

We obtained Figure 5 by considering the spinodal stability of our two component system. 
The spinodal points are found from the following equation 



J{<i>As) = 



= (36) 



with T — F{T^q)v/ {kTV) being the reduced equilibrium free energy. In Figure 5 in the area 
above the curve for a certain x the equation J(0, 0^ = 0.1) = has two solutions, while in 
the area below the curve it has no solutions for physical values of (j). It should be noted that 
in our two component system this condition on the existence of spinodal at a given salt is 
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strictly speaking not equivalent to existence of phase separation (to demand that the critical 
point be at (pg — 0.1 is yet another different condition). However, in the considered system 
dissolution occurs at very low 4>, so the three conditions yield numerically close results. 

IV. CONCLUSIONS 

Complexation in solutions of oppositely charged polyelectrolytes can be accompanied by 
thermoreversible crosslinking of oppositely charged monomeric ions. A local electrostatic 
binding energy can exist between oppositely charged units when they tend to be dehydrated 
in the vicinity of each other or due to non-classical specific interactions. In this work we 
have investigated how the three different mechanisms (long-range electrostatics, crosslinking 
and backbone hydrophobicity) define the properties of polyelectrolyte complexes at different 
salt concentrations. 

In our approach we obtain self-consistently the fraction of crosslinked charged monomers 
(conversion) and the Debye-Hiickel collective fluctuations contribution to the free energy 
(which depends on conversion). Accordingly, the degree of conversion is determined both 
by the local binding energy and by long-range electrostatics. We find that the long-range 
charge fluctuations always promote crosslinking. Given that the magnitude of the Debye- 
Hiickel contribution decreases with increasing salt, the fraction of crosslinked monomers also 
monotonically decreases with increasing salt. 

The polymer concentration in the precipitate is largest at low salt concentration, when 
the screening of interactions between monomeric ions is weakest. The complex concentration 
generally decreases monotonically with increasing salt concentration. The rate of complex 
dilution with addition of salt and the concentration of monomers in the precipitate at high 
salt are strongly dependent on the value of the van der Waals attractions and on the binding 
energy. Non-selective net van der Waals attraction between the monomers of both the 
positively and negatively charged chains enhances the complexation in a way broadly similar 
to crosslinking due to local binding energy. The dilution is very rapid and the monomer 
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concentration of the complex goes to zero in the case of zero binding energy and zero net 
van der Waals attraction (good solvent condition) . Instead, for large values of either type of 
the short-range attractions, as the salt concentration increases the monomer concentration 
in the complex generally nearly saturates to a nonzero constant value (or slightly increases 
for sufficiently large x)- Our results are consistent with experimental observations^^"^^ in 
which, depending on the type of polymers used, with the addition of salt the complexes 
can either dissolve, or their density can remain stable. In some cases initial dissolution 
and subsequent re-entrant precipitation is observed^^, the type of behavior obtained in our 
theory for marginal solvent. 

There is an important competition between complexation due to charge fluctuations 
and non-linear thermoreversible linking, which is especially interesting in good solvent con- 
ditions. Unassociated charged monomeric groups induce complexation due to long-range 
electrostatics. However, once they form crosslinks they practically do not contribute to 
long-range attraction. Therefore at high conversion rates complexation is mostly due to 
effective crosslinking attraction in an effectively neutral polymer solution. This competition 
leads to an interesting non-monotonic behavior of monomer concentration in the precipi- 
tate with increasing the non-linear binding energy in the case of zero salt concentration: 
with the increasing binding energy the monomer concentration in the precipitate passes 
through a minimum. Remarkably, the variation in the complex density is rather small, 
that is, crosslinking and long-range electrostatic attractions give rise to complexes of similar 
density. 

Another important effect in polyelectrolyte complexation is the difference in salt con- 
centration inside and outside the precipitate. When the non-linear binding energy is small 
the difference in salt concentration in and out of the precipitate is negligible. However, for 
large binding energies (or large values of x) this difference rapidly grows as the overall salt 
concentration increases, the precipitate being depleted of salt due to increasing importance 
of hardcore interactions.^^ We find that for large binding energies and/or in a bad solvent 
the difference between salt concentrations in the complex and in the bulk has a significant 
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effect on tfie density of the complex at high salt concentrations. An interesting limit to an- 
alyze includes the addition of non-linear correlations among the ion pairs when the fraction 
of charged units increases as in the case of strongly charged chains in oppositely charged 
multivalent ion solution^^'^^ where even denser precipitates are expected. 

Acknowledgments. This work was supported by NSF grant DMR 041446. A. E. 
acknowledges partial financial support of NSF grant EECO 118025. 

APPENDIX A: CALCULATION OF THE STRUCTURAL CORRELATION 

FUNCTIONS 

In this section we calculate the structural correlation functions of the solution of associ- 
ating oppositely charged polyelectrolytes. Charged groups on chains are treated as stickers 
that can associate with ionic groups of an opposite sign. We consider only the symmetric 
case of oppositely charged homopolymers of equal degree of polymerization N — Ni — N2, 
present in solution with equal concentrations pi — p2 — p/2. (Indices 1 and 2 refer to 
positively and negatively charged chains, respectively.) We also assume that the chains have 
the same chemical structure, that is the same bond length b = bi = 62, and the distance 
between charges a — bf~^^^ (/ is the fraction of charged monomers on the chains). We 
assume Gaussian statistics for chains. 

The expression for the structural correlation function (9) between two different types of 
monomers a and /3 reads (for detailed derivation see ref 23) 

^a/3(g) = Ep(C)^S3(?) (Ai) 

C 

Summation in (Al) is over all topologically different clusters formed due to association of 
polymers, with p{C) being the number concentration of a cluster having structure C and 
Qa^iq) the molecular structural correlation function. In (A2) the summation runs over all 
monomers of types a and (3 of the cluster C. The average is over the conformations of the 
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cluster, it can be written as 



- ^""'IfVi^^'^" = V /e-^"'-^'/a(rc)rfr, (A3) 



where /^(rc) is the probabihty function of finding the cluster in conformation F^, and the 
integration is over the entire configurational space of the cluster C. 

To calculate the correlation functions (Al) we employ the grand canonical diagrammatic 
technique described in ref 23. (For details of the diagrammatic technique as well as appli- 
cations to other related systems see refs 22, 23, 47.) As is shown in ref 23 the correlation 
function can be expressed as the sum of all two-root diagrams 

where z is the fugacity of the chain {z = exp(— /^//cT), /j, is the chemical potential), W{C^^) 
is the statistical weight of the cluster C^^ with two marked monomers of types a and /3, 
and S{C^^) is its symmetry index. 

In order to calculate Qapio) it is convenient to introduce the following generating func- 
tions. Let us introduce the generating function of all one-root diagrams t (diagrams with 
one marked monomer), which can be calculated recursively as 

t^l+uzt''-^^ (A5) 

Here to is the statistical weight of the crosslink: to — exp(£), where e is the absolute value 
of the dimensionless crosslink bond energy £ = jkT . Now we can write the sum of all 
labeled diagrams with two labels belonging to the same chain as 

where we have introduced the correlation function of one homopolymer chain g'(g). Note 
that the two labeled points are free, that is they have no diagrams attached to them. We 
need to introduce also a closely related to Eg sum of all diagrams where the two labels 
cannot belong to the same monomer, 
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/•(<;) = 9(9) - 1 (A8) 



Since the system is symmetric we need to calculate only two functions gu — g22 and 
912 — 921- Using definitions (A5-A8) the autocorrelation function 9u{(l) can be written as 
the following series 

9n{q) = t^^g + f^giu'JlhWJ:^ + ^^^^(a;'^;,) VSg + ... (A9) 

Here, the first term is the sum of all diagrams in which the two labels belong to the same 
chain. Next terms result from summation of all diagrams with labels belonging to different 
chains. In our model only oppositely charged chains can associate with each other. The 
second term is the sum of all diagrams in which the two labels (marking monomers on chains 
of type 1) are separated by a chain of opposite charge (type 2). The next term in (A9) comes 
from summation of all diagrams with the insert between the labeled chain comprised of three 
chains (sequence 2 — 1 — 2). Higher terms correspond to summation of diagrams with a higher 
number of chains in the insert between the labeled chains. 

In (A9) we introduced the effective statistical weight of the crosshnk 

co' = coe-i'"'/'' (AlO) 

which takes into account the correlations of monomers in a crosslink, which we assume to be 
Gaussian (b is the bond length). Note that since one monomer can form only one crosslink 
the presence in a diagram of a chain connecting the two labeled chains corresponds in eq 
A9 to a generating function E/j , in which summation runs over i ^ j. It is easy to see that 
the term in (A9) corresponding to the sum of all diagrams separated by 2n — 1 chains reads 
i^Eg(a;'E/i)^"~^a;'Eg. It is easy to sum the infinite series (A9) as 



l+U;'T,gLj'{T,g - Eft) 

1 - (u^'^nf 



9iM = . \ )J.2 ' (All) 

Analogously for 9i2{q) we obtain 
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gi2{q) = t'^J^gCu'j^g + eUgiu'Uhfu'llg + ^11 g{u' h)^ U' g + ... (Al2) 

where the first term corresponds to all diagram with the labels belonging to two neighboring 
chains (of opposite charge), the next term — all diagrams with two chains separating the 
labeled chains and so on. The series (A12) can be summed as 

9i2{q)=eT.g ^^^i (A13) 

To use the correlators (All) and (A13) in the free energy we need to change from the 
variables z and t to the number concentration of charged monomers pi (for our case pi = P2) 
and conversion F. Conversion F is defined as the fraction of charged monomers in crosslinks 

where p^^^ is the number of crosslinks. As is shown in refs 22, 23 the concentrations are 
given by 

zt^ 

Pi = -i^N (A15) 
p(^)=a;Q^Art^-^y (A16) 

Combining it with (A5) and (A6) we obtain 

t%=pig(g) (A17) 

uj'T.h = Te-'''''''%{k) = T'h{k) (A18) 

which substituted into (All) and (A13) yields 

= Pigiq) ^'^^"^^ o (A20) 

Note that in (A18) we introduced the effective conversion F', which takes into account the 
correlations of ions in the crosslink. The divergence of the correlators at T'h{q = 0) = 1, i.e. 
F(A^ — 1) = 1 corresponds to gelation. 
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The value of F can be obtained from the definition of F (given by eq A14) using relations 
(A5), (A15), and (A16). Conversion turns out to be determined by the unmodified mass 
action law 

jY^Yf " ^'"^ " ^'""^ ^^^^^ 

This reflects the fact that in this Appendix we considered an ideal associating system (no 
other interactions except association are present). Conversion for the solution of associating 
polyelectrolytes is obtained from the minimization of the total free energy, which includes 
interaction terms. The resulting modifled mass action law (28) differs from eq A21 in that 
it has a long-range electrostatic contribution to the effective binding energy. 
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Figure Captions 

Figure 1. (a) Polymer volume fraction in the precipitate as a function of the salt 
volume fraction in the supernatant (pg for different binding energies e. Dash lines correspond 
to the assumption of equality of salt concentrations in the precipitate and supernatant: 
(p^f^ = (ps- (b) Conversion T for curves of plot (a), (c) The difference of salt volume fractions 
in the precipitate and supernatant 0^^^ — (ps for the curves of plot (a). 

Figure 2. Effect of the Flory-Huggins x-parameter on the change of polymer volume 
fraction in the precipitate with increasing salt in the system (pg. 

Figure 3. Variation of the density of the precipitate for salt- free system with changing 
binding energy e. Different curves correspond to different values of the x-parameter. 

Figure 4. Coexistence lines for phases with different polymer volume fractions at a 
given salt volume fraction (pg- The effect of varying chain length N is illustrated. 

Figure 5. Stabihty of the precipitate to addition of salt. For all values of the bond 
energy e and the chain length N below the curve for a corresponding value the precipitate 
dissolves when the salt concentration is increased beyond 4>s > 0.1. 
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